Build up of off-diagonal long-range order in microcavity exciton-polaritons across the 

parametric threshold 
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We report an experimental study of the spontaneous spatial and temporal coherence of polariton 
condensates generated in the optical parametric oscillator configuration, below and at the parametric 
threshold, and as a function of condensate area. Above the threshold we obtain very long coherence 
times (up to 3 ns) and a spatial coherence extending over the entire condensate (40 /mi). The very 
long coherence time and its dependence on condensate area and pump power reflect the suppression 
of polar iton-polariton interactions by an effect equivalent to motional narrowing. 
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Planar microcavity exciton-polaritons are composite 
particles generated by strong coupling between excitons 
and cavity photons that follow the Bose statistics. They 
have been extensively investigated in the past decades 
for their unique properties, such as very light effective 
mass and the consequent high critical temperature for 
condensation pQ. It is well known that the occurrence 
of such a phase transition is accompanied by the onset 
of long range phase coherence [2]. In fact, since the un- 
ambiguous demonstration of polariton condensation its 
coherence properties have been investigated [3]. More- 
over, the extension and eventually the decay of the spa- 
tial coherence can give useful information on the kind of 
transition the system undergoes [4j [5] . Recently, polari- 
ton condensates have attracted intense attention for their 
potential application in the field of quantum information, 
and also some properties suitable for such applications, 
such as Josephson oscillations and spin-switching, have 
been already demonstrated [6[ [7] . Therefore the achieve- 
ment of extended spatial and temporal coherence for the 
polariton condensate is crucial. 

Here we study the spatial and temporal coherence 
properties of nonequilibrium polariton condensates un- 
der an optical parametric oscillator (OPO) regime. The 
OPO is a third order non-linear phenomenon, arising 
from interactions between the excitonic components of 
the polaritons. It takes place when two pump polari- 
tons at the inflection point of the lower polariton branch 
(LPB), scatter efficiently into a signal and an idler po- 
lariton [8j |9]. Spatial coherence properties have been 
studied theoretically for the OPO condensate by Caru- 
sotto and Ciuti [10]. They investigated numerically the 
first order coherence function (g^) for a finite conden- 
sate paying special attention to its behaviour across the 
OPO parametric threshold {Eth) [ZD- H is found that, 
for excitation frequencies uj p below that of the threshold, 
has a finite correlation length. Increasing uj p in or- 
der to approach Eth-, mantaining fixed the pump angle 



(and therefore the wavevector), they predict the build 
up of macroscopic phase coherence extending over the 
entire condensate [ID] . In this regime the spatial fluctu- 
ations are negligible, so the temporal coherence proper- 
ties should be captured by the theory of Whittaker and 
Eastham [12]. In this theory the temporal coherence is 
limited by fluctuations in the particle number, which due 
to the polariton-polariton interactions imply a broaden- 
ing of the emission [13-15 . This broadening mechanism, 
however, could be suppressed by motional narrowing [12], 
so that for appropriate pump powers and condensate ar- 
eas very long coherence times could be obtained. De- 
spite the fact that long coherence times and extended 
spatial coherence have been predicted theoretically for 
the OPO configuration, only now our experiments con- 
firm such predictions. 

One of the main factors that limits the coherence is the 
quality of the cavity: the presence of defects creates a 
disorder potential that traps the condensate, potentially 
leading to multi-mode and inhomogeneous states [TP] . 
Another, very important, detrimental effect is caused by 
the fluctuations of the excitation laser that hinder the 
attainment of the intrinsic coherence of the condensate. 
Many studies [3j [16] [17] have been performed by non- 
resonantly pumping the microcavity, in such cases the 
resulting distribution of the population at the bottom of 
the lower polariton branch is subjected to fluctuations 
due to the reservoir of particles at the bottleneck. These 
fluctuations, broadening the distribution of polaritons in 
energy and momentum space, translate, according to the 
Wiener-Khinchin identity [I8j [19] , into a faster decay of 
the temporal and spatial coherence. Although analogous 
broadening by the fluctuating population of pump polari- 
tons can occur in the OPO [15], the threshold pump den- 
sity (Pth) for a resonant-gain process is much lower than 
that for non-resonant gain and the effect of the reser- 
voir polaritons is less important. Thus the coherence ex- 
hibited by an OPO polariton-condensate is expected to 
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FIG. 1: (color online) (a) Typical real space emission 
intensity (color scale) from the signal condensate 
generated by parametric scattering, at detuning 
S ~ —3.8 meV and pump power P = 10kW/cm 2 . (b) 
Corresponding energy spectrum of the signal emission 
(spectrometer limited) . 



decay over much longer times and larger distances than 
the one produced by non-resonant techniques. To avoid 
that the pump fluctuations conceal the genuine coherence 
properties of the condensate we use a CW monomode 
laser with a very narrow bandwidth of 75 kHz to excite 
the sample, which is a high-quality A-microcavity grown 
by molecular beam epitaxy, composed by 16 periods of 
A/4 AlAs/Alo.isGao.ssAs on its top and 25 periods on 
its bottom. A single 10 nm GaAs/Alo.3Gao.7As QW is 
placed at the antinode of the cavity. The Rabi splitting 
is 4.2 meV and we select a detuning of 5 ~ —3.8 meV. 
In order to investigate the spatial coherence properties 
of a true 2D condensate we create one with a diameter 
of ~ 40 /im, by pumping resonantly the LPB close to its 
inflection point, with a Gaussian laser beam, maintain- 
ing the temperature of the sample at 10 K. We adopt 
the experimental conditions to drive the sample into the 
OPO regime at the parametric threshold, Eth, an d then 
we follow a similar approach to that in [ID] , i.e., we de- 
crease the pump frequency uj p keeping fixed the angle of 
incidence. Typical results for the OPO signal are shown 
in Fig. [I] They reveal a bright and rather homogeneous 
emission in real space (Fig. []Ja)), with a spectrometer- 
limited emission spectrum (Fig. [TJb)). 

Therefore, to evaluate the coherence of the conden- 
sate, since its true energy linewidth is masked by our 
spectrometer, the signal emission is analyzed with the 
help of a Mach-Zehnder interferometer. On one of the 
arms of the interferometer a retro reflector mirror is 
mounted, which flips the image of the condensate in a 
centrosymmetric way. Recombining, at the output of 
the interferometer, the image of the condensate with its 
flipped image allows evaluating the coherence of each 
point (x,y) with respect to its opposite counterpart 
(— x, —y). In this way we access to the first order coher- 
ence function gW [(x,y) , (-x, -y)\ = ( ^ ( i^jfi^_X-y)) 
[2j [3]. To extract the coherence of the condensate the 
retroreflector mirror is moved, using a piezo translation 



stage, by a few wavelengths around the zero delay posi- 
tion, adding a short delay At to the retroreflector arm 
with respect to the other arm. For each measurement 
the intensities of the two arms were acquired and used 
in the following normalized expression of the interfer- 
ence pattern: I Nor m [(x,y) , (-x, -y)] = ^oi-Irr-Im = 

^ V 1 RR 1 M 

[(#, y) , (— x, —y)} sin (At + cpo) where Irot is the to- 
tal intensity of the interference pattern, Irr (Im) is the 
intensity recorded from the arm with the retroreflector 
(single mirror), and cpo is the initial phase of the con- 
densate. For each point of the condensate a sinusoidal 
evolution of the normalized interference as a function of 
the delay At is recorded. The first order coherence of 
each point of the condensate is extracted with the help 
of the previous formula; and repeating this procedure for 
each point we can reconstruct the entire coherence map. 

At EtHi the interference pattern extends along the en- 
tire condensate (Fig. |2ja)) and the corresponding co- 
herence map (Fig. (2^b)) reveals that a constant phase 
is developed through its entire area, validating the the- 
oretical predictions [ID]. This is shown in detail in Fig. 
|2|c), which depicts the coherence of a horizontal profile 
(0.26 /im width) taken at the center of the map. The fact 
that the degree of coherence is almost flat along the en- 
tire condensate demonstrates that only a single coherent 
mode is formed, in contrast with what has been found 
in previous works [TDJ [20] where disorder of the sample 
caused the development of several modes. The extension 
of the coherence here is ~ 40 /im, the largest reported 
until now in 2D microcavity polaritons, to the best of 
our knowledge; only in ID microcavities, for propagating 
condensates, a larger spatial coherence has been shown 
[2l] , It is worth to note that in an OPO process, the co- 
herence properties of the signal are not inherited from the 
excitation laser [22]. Although the theory of an infinite 
2D condensate predicts the absence of long-range order, 
and a power-law decay of the spatial coherence [4], in a 
finite dimension system such as ours, determined by the 
pump laser size, a constant coherence of the condensate 
can be obtained [3 . 

By decreasing the pump energy by SE = 0.01 meV to- 
wards lower energies, leaving fixed the pump angle, the 
phase matching conditions for the OPO are not fulfilled 
anymore, although an emission is obtained at the bot- 
tom of the LPB. The interferometric analysis of the emis- 
sion from the cloud of non-condensed polaritons gives an 
interference pattern only at the center of the emission 
(Fig[2jd)). The corresponding coherence map (Fig[2je)) 
reveals a ~ 7 /im coherence length, of a few polariton de 
Broglie wavelengths, extracted by fitting the profile of the 
spatial decay (Figj2jf)). These results demonstrate the 
predictions of Ref. [ID] about the coherence of polariton 
emission below and at the condensation threshold. 

Measuring g^ at a given point for different time delays 
r permits the evaluation of the signal temporal coherence 
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FIG. 2: (color online) Interference pattern (a/d) and corresponding coherence (b/e) of the condensate generated 
at/below the parametric threshold, Ethi corresponding to a pump energy E\ = 1552.57 meV/ E2 = 1552.56 meV 
and power density P — 15.8kW/cm 2 . Horizontal profiles at the center (c/f) showing the constant coherence along 
the entire condensate for E\, and an exponentially decaying one for E2. 



gi 1 ) ( r ) = ^i^i^ll) . i n our case we obtain for dif- 
ferent values of the condensate area A and pump power 
density P, as shown in Fig. [3j We obtain a predomi- 
nantly exponential decay for (r), with a coherence 
time reaching T c = (3.2 ± 0.8) ns. This is much longer 
than the largest value reported in the literature until now 
for the same material system [T5] . 

Despite the lack of spatial fluctuations, the coherence 
time of our condensate, though long, is clearly finite. In 
general, such a finite correlation time in a state with per- 
fect spatial order is caused by finite-size fluctuations [23] , 
and reflects the absence of true phase transitions in finite 
systems. A well-known example is the Schawlow-Townes 
form for the coherence time of a single laser mode with an 
average of N photons, T c oc N. Since N oc A when the 
control parameter P is constant, the Schawlow-Townes 
result implies the scaling form T c (x A. However, as 
shown in detail below, this scaling law is violated by more 
than an order of magnitude for our condensate, and our 
results cannot be understood as a straightforward effect 
of increasing the particle number. Our results may be in- 
terpreted in terms of the theoretical model of Ref. [12] , in 
which the linewidth arises from fluctuations in the num- 
ber of particles. Due to the interactions, such number 
fluctuations imply energetic fluctuations of the emission, 
leading to a broadened line or a decay of g^ (r). This is 
captured by Kubo's result 
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for the emission from a transition whose energy fluctu- 
ates, where r c is determined by the width of the (generi- 
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FIG. 3: (color online) (a)Temporal coherence decay 
above parametric threshold for condensate areas 
A = 70/im 2 (blue dots) and 2.7 x 10 3 /im 2 (black 
diamonds) at a pump power density P = 10kW/cm 2 ; 
(b) Temporal coherence decay for P = 5.4kW/cm 2 
(black squares) and 10kW/cm 2 (blue dots) at 
A = 70/im 2 . Lines are fits to Eq.(l). 



cally Gaussian) distribution of the energy level and r r is 
the characteristic timescale on which the fluctuations oc- 
cur. If r r ^> r c the fluctuations are slow and the emission 
reflects the energy distribution of the level, so that the 
coherence decays on timescale r c . If however r r <C r c a 
motional narrowing effect leads to an exponential decay, 
with a much longer coherence time T c = r 2 /2r r . 

To elucidate the effects of number fluctuations, we fit 
the data in Fig. [3] to the Kubo formula 0. Pumping 
with the same power density P = 10kW/cm 2 , the con- 
densate with a larger area A = 2.7 x 10 3 /im 2 presents 
r c = (0.97 ±0.11) ns and r r = (0.15 ± 0.02) ns fulfilling 
the inequality for the motional narrowing regime and giv- 
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ing a coherence time T c = (3.2 ± 0.8) ns. With a smaller 
area A = 70/mi 2 , we find a decrease of both parameters 
t c = (0.20 ± 0.08) ns and r r = (0.023 ± 0.007) ns, still 
indicating the system being in the motional narrowing 
regime but giving a shorter T c = (0.9 ± 0.4) ns. As well 
as moving to larger areas, another way to have longer 
decay times suggested by [12 is to increase the pump 
power. The results shown in Fig. [3^b), corresponding to 
two different powers keeping now the size of the conden- 
sate constant, also demonstrate that increasing the power 
a longer coherence time is achieved. In fact the coherence 
time for a condensate of relatively small area A = 70 fim 2 
goes from T c = (0.5 ± 0.3) ns at P = 5.4kW/cm 2 , to 
T c = (0.9 ± 0.4) ns at P = 10kW/cm 2 . We note that 
the coherence decay of the smallest condensate at lowest 
power, black squares in Fig. [3^b), has considerable struc- 
ture and does not follow the Kubo form, suggesting that 
this condensate is not a single mode. From the Kubo 
fits to our data, for a given power, we can extract the 
scaling r c cx A 0A±0A . This is consistent with the domi- 
nant dephasing mechanism being the polariton-polariton 
interaction, suppressed by motional narrowing, in which 
case the Kubo formula should hold with r c cx A - 5 [T2] . 
Although this implies an increase in the coherence time, 
T c , with area, our results show that this increase is re- 
strained because r r also increases, so that motional nar- 
rowing becomes less effective. 

The relatively small size of our system and the use 
of pump powers which are close to threshold, Ptk, are 
responsible for the observed increase of r r with area. 
In the thermodynamic limit, A — >> oo, the occupation, 
n, of a single mode with linear gain 7, linear loss j c 
and nonlinear gain T obeys the mean-field rate equation 
dn/dt = (7 — 7 C — Tn)n [24]. Linearizing we see that the 
damping time for number fluctuations is r r = I/I7 — r y c \. 
This timescale is therefore independent of area when 
the rate equation is valid. However, at power thresh- 
old 7 = 7c and the rate equation predicts a divergence 
in the relaxation time, which in the finite system must 
be cut off by fluctuations. Thus, in the threshold re- 
gion, r r initially grows with area, as we observe, be- 
fore eventually saturating. We obtain a similar behav- 
ior from numerical solutions to the dynamical model de- 
scribed in Ref. [12] in the threshold region. If we as- 
sume a polariton lifetime tq = 3 ps, then our results for 
P = 10kW/cm 2 , A = 70 fim 2 imply the critical slowing 
factor r r /ro ~ 8. Fitting this to the dynamical model at 
threshold gives a resonable value for the gain saturation 
density, p s = n s /A — 230/(70/im 2 ). Since p s is inde- 
pendent of A, we use this value in the model to predict 
r r = 140 ps for our largest condensate, consistent with 
experiments. 

For completeness we also studied the temporal coher- 
ence (not shown) under two further conditions: a) Excit- 
ing the LPB at the inflection point, for power densities 
slightly above but very close to Pth-, we find a Gaus- 



sian decay of g^\r), as previously observed for conden- 
sates with non-resonant pumping [16] . This is consistent 
with the expected increase in r r approaching Pth caus- 
ing a crossover from the motional-narrowing to the static 
regime, b) Exciting far from the inflection point (no con- 
densation achieved) we observe a very fast exponential 
decay of coherence on the range of the polariton lifetime. 

In summary, we have investigated the spatial and 
temporal coherence properties of polariton condensates 
generated by parametric scattering. The spatial coher- 
ence was studied below and above the OPO condensa- 
tion threshold, obtaining the largest value for the coher- 
ence length measured for 2D GaAs microcavity polari- 
tons. Measurements of the temporal coherence reveal 
a predominantly exponential decay caused by polariton- 
polariton interactions in a motional narrowing regime. 
Although similar to the exponential decay associated 
with the Schawlow-Townes linewidth, the mechanism is 
completely different, and gives a different dependence of 
coherence time on condensate area. By constructing con- 
densates of large area we are able to achieve long coher- 
ence times, which is a crucial step forward for exploiting 
polariton condensates in quantum and ultrafast devices. 
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